Building micro-soccer-balls with evaporating colloidal fakir drops 
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Evaporation- driven particle self-assembly can be used to generate three-dimensional microstruc- 
tures. We present a new method to create these colloidal microstructures, in which we can control 
the amount of particles and their packing fraction. To this end, we evaporate colloidal dispersion 
droplets on a special type of super hydrophobic micro- structured surface, on which the droplet re- 
mains in Cassie-Baxter state during the entire evaporative process. The remainders of the droplet 
consist of a massive spherical cluster of the microspheres, with diameters ranging from a few tens 
up to several hundreds of microns. We present scaling arguments to show how the final particle 
packing fraction of these balls depends on the dynamics of the droplet evaporation. 
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Evaporation-driven particle self-assembly is an ideal 
mechanism for constructing micro- and nanostructures 
at scales where direct manipulation is impossible. For 
example, in colloidal dispersion droplets with pinned con- 
tact lines, evaporation gives rise to the so-called coffee- 
stain effect [ ]: a capillary flow drags the particles to- 
wards the contact line to form a ring-shaped stain. Such 
a flow not only aggregates the particles, but is also able to 
organize them in crystalline phases [ - ]. Similar mech- 
anisms such as the "convective assembly" [ , ] are cur- 
rently successfully used to produce two-dimensional col- 
loidal crystal films. To obtain three-dimensional clusters 
of micro-particles, colloidal dispersion droplets which are 
suspended in emulsions [ -10] or kept in Leidenfrost levi- 
tation [ ] are used. With these techniques, new colloidal 
structures arise from the geometrical constraints during 
the drying [ ]. This problem of organization of parti- 
cles into a spherical topography dates back to the days 
of the first models of the atom and has been extensively 
studied by Bausch et al. [13, 14]. The main drawback of 
these three-dimensional assembly techniques, however, is 
the lack of control on both the amount of particles and 
the particle arrangement in the remaining structures. 

In this work, we devise a new, controlled way of 
generating on-demand self- assembled spherical micro- 
structures via droplet evaporation on a superhydropho- 
bic surface (see figure 1). We present scaling arguments 
to predict the particle arrangement in the microstruc- 
tures formed, based on the dynamics of the evapora- 
tion process. To generate the microstructures, we evap- 
orate colloidal dispersion droplets on a special type of 
superhydrophobic substrates. In most of the cases, a liq- 
uid Cassie-Baxter state drop evaporating on a superhy- 
drophobic surface will eventually suffer a wetting tran- 
sition into a Wenzel state, i.e. it will get impaled into 
the structure and loose its spherical shape [15, 16]. Here, 
however, we use a surface that combines overhanging pil- 
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FIG. 1. (a) A droplet of colloidal solution is left to evaporate 
on a superhydrophobic surface. As the solvent evaporates, the 
particle packing fraction increases. Once all the solvent has 
completely evaporated the colloidal particles have aggregated 
to form a spherical particle conglomerate: a colloidal supra- 
ball. (b) Top view of the resultant compact colloidal supraball 
left on the superhydrophobic surface after evaporation. The 
micropillars forming the structure are seen as circular objects 
around the supraball. 



lared structures [ ] with a hierarchical nano-structure 
(figure 2c). These surface properties impose a huge en- 
ergy barrier for the wetting transition to occur, and 
therefore the droplet will remain almost floating over the 
structure in a Cassie-Baxter state during its entire life 
[18]. A typical result can be observed in figure 1: a water 
droplet containing ljam polystyrene particles (concentra- 



tion 0.08% weight and initial volume 5/iZ) evaporates on 
the superhydrophobic surface at room temperature and 
30% humidity. After a typical evaporation time of 45 
minutes, the solvent is completely evaporated and only 
the colloids are left upon the substrate. Remarkably, 
the particles (~ 10 T in this particular case) are not just 
lying scattered over the substrate but they have aggre- 
gated and form a spherical macro-cluster resting on top 
of the micro-pillars, which we call colloidal supraball. We 
do not observe shell formation and buckling during the 
evaporation of the droplets [ ]: the supraballs we ob- 
tain are solid, and present a high mechanical resistance 
and stability. When looking closer at the surface of these 
particular colloidal supraballs as shown in figure 2d, one 
can identify crystalline flat patches which resemble the 
pentagonal patches present in a soccer ball. 

To understand the final structure of these supraballs, it 
will turn out crucial to understand the dynamics of the 
droplet evaporation. The fact that we do not observe 
shell formation, suggests that the particles do not influ- 
ence the droplet evaporation. To test this hypothesis, 
we compare the evaporation dynamics to that of a liquid 
drop that does not contain any particles. The evapora- 
tive mass loss from such a drop is typically governed by 
the diffusion of vapor molecules in the surrounding air 
[1, 19, 20]. For diffusion-limited evaporation, the rate of 
volume change of the drop is given by 



dV 
~dt 



D'K 



(1) 



where R is the drop radius, and D' = D va Ac/p, with 
D va the diffusion constant for vapor in air, Ac the vapor 
concentration difference between drop surface and the 
surroundings and p the liquid density [ ] . One might have 
expected the evaporation rate from the drop surface to be 
proportional to the droplet surface area ~ R 2 . However, 
the vapor concentration gradient is proportional to l/i2, 
and therefore the total evaporation rate is proportional to 
R [21]. If the droplet evaporates with a constant contact 
angle, we find that, since V ~ R 3 , 



the interface speed increases dramatically as the droplet 
shrinks and the maximum speed reached in the exper- 
iment will be determined by the final radius Rbaii- As 
we will show further on, this increase in interface speed 
determines the particle packing inside the supraballs. 




FIG. 2. (a) Tilted view of the supraball in contact with the 
microstructure. (b) Detail of the contact area, (c) Magni- 
fied view of the micropillars forming the microstructure. (d) 
Close-up of the supraball surface. The distribution of crys- 
talline patches resemble the pentagons in a soccer-ball. 

The final size of the ball depends on the number of par- 
ticles inside the drop. This can be tuned by manipulat- 
ing either the initial particle concentration or the droplet 
size. In our experiment, the ball size was in the range 
100 < Rbaii/Rp < 1000, with R p the particle radius. 
Clearly, the exact final size of the ball will not only de- 
pend on the amount of particles in the system but also 
on their packing fraction. We define the packing fraction 



R{t)~[D'{t f -t)\ 



1/2 



(2) 



Here tf is the total droplet lifetime in case no particles 
are present, for which the drop radius reaches zero. In the 
present case the drop radius saturates at a finite radius, 
Rbaih at a time i = tf — R^ all /D' , corresponding to the 
moment where the particles become densely packed. In 
figure 3 we plot the colloidal droplet radius versus tf — t. 
Our experimental data for different number of particles 
are in very good agreement with the 1/2-power law. This 
confirms that the particles do not influence the evapora- 
tion process, until the final radius Rbaii is reached. The 
scaling (1) implies that the speed with which the interface 
is moving inwards, is given by dR/dt ~ D' /R. Hence, 



*_/v< Rp 



Rball 



(3) 



where TV is the total number of particles in the droplet. 
The final supraball radius Rbaii is accurately determined 
from Scanning Electron Microscope (SEM) images. If 
the packing fractions were identical for all supraballs, one 
would expect that Rb a ii/R P ~ A/ -1 / 3 , as depicted in figure 
4a [ ]. However, in figure 4b we see that the packing 
fraction strongly depends on the number of particles in 
the system. 

As the number of particles increases, the packing frac- 
tion approaches that of a perfect Hexagonal Close Pack- 
ing configuration, in which case one would find $ = 0.74 




(tf-t)(s) 



FIG. 3. The radius of the droplet plotted against tf—t with tf 
the lifetime of the droplet and t the actual time. The triangle 
indicates the 1/2-power law, the dots represent the data sets 
for 7 different experiments, where the number of particles was 
varied. For a certain t < tf the final ball size Rbaii is reached. 
The final time was extrapolated as tf — i + R% a u/D'. 
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FIG. 4. (a) Supraball to microparticle diameter Rb a ii/R P 
plotted against the total amount of particles N in the sys- 
tem, (b) Packing fraction $ strongly depends on N. Blue 
dots represent experimental measurements and the red solid 
line corresponds the most efficient particle packing <E> = 0.74 
(hexagonal close packed), the dashed line marks $ = 0.64, 
random close packed. N c is the critical number of particles, 
above which we find an ordered ball structure 



[23], hence, we have an ordered particle packing inside the 
balls. On the other hand, the supraballs with a smaller 
amount of particles show remarkably low packing frac- 
tions, even below the Random Close Packing (RCP) limit 
(<1> = 0.64) [ ], corresponding to a disordered particle ar- 
rangement. The balls which show packing fractions below 
the RCP limit contain several empty cavities. Remark- 
ably, the final configuration reached, seems to depend on 
the number of particles in the system. In figure 4b we 
indicated the critical number of particles N c « 3 • 10 6 
when the packing fraction reaches that of a RCP. For 
N < N c we get a loose, disordered particle packing in the 
supraball, whereas for N > N c we get a densely packed, 
ordered supraball. 

What causes the transition from ordered to disordered 
packings, and what determines the critical number of 
particles? To explain this, we follow a similar approach 
as in Marin et al. [ ]: we compare the time scale on which 
particles can arrange by diffusion to the hydrodynamic 
time scale for the particle transport by convection, given 
by the inward motion of the liquid-air interface. If the 
diffusion time is small compared to the hydrodynamic 
time, particles can arrange into an ordered packing. The 
diffusive time scale is td = Rl/D, with R p the particle ra- 
dius and D the diffusivity of the particles in the liquid [ ] . 
The hydrodynamic time-scale is t^ = Lj | ^ | . Here R(t) 
is the droplet radius, and L is the typical inter-particle 
distance, which depends on the particle concentration as 
L = N~ l ^R, as long as the solution is dilute (L ^> R p ). 
We define the ratio of both time scales as: 



A(t) 



t h 



dR(t) 



dt 



td 
L 



D \R(t) 



(4) 



where in the last step we used (1) to replace dR/dt ~ 
D'/R. 

From (4) we observe that A(t) increases as the droplet 
radius becomes smaller during the evaporation (see figure 



3), until the limit R = Rball is reached. A cross-over be- 
tween the time-scales is reached when the hydrodynamic 
time becomes equal to the diffusion time, hence when 
A = 1. If the cross-over is reached when R ^> Rbaiu the 
amount of crystalline clusters is still very small. From 
this point in time onwards the interface speed is too high 
for the particles to further arrange in a crystalline way 
[ ]. Instead, they are pushed together in a random ar- 
rangement, with a low packing fraction. If the cross-over 
is reached when R < Rbaiu the particle packing is already 
dense and ordered, and we find a high packing fraction. 
For all droplets the evaporative mass loss, and hence the 
decrease in radius, is the same (see figure 3), hence, the 
moment when the particle packing becomes sufficiently 
dense for particles to arrange depends solely on the num- 
ber of particles in the droplet. If N is high (N > iV c ), 
this moment is reached relatively early, i.e. well before 
.4 = 1, and we get an ordered particle packing inside the 
supraballs. Using that Rball/ Rp ~ N 1 / 3 and considering 
A = 1, we find from (4) the critical number of particles 
above which we obtain ordered supraballs 



N r . 



(5) 



This result emphasizes that the transition is governed 
by two diffusion processes: the diffusion of vapor, deter- 
mining the speed of evaporation, versus the diffusion of 
particles inside the drop. The ratio of diffusion constants 
selects the critical number of particles. In our exper- 
iment D' = 3 x lO" 10 m 2 /s and D = 2 x lO -13 m 2 /s, 
from which we find that 7V C ~ 10 9 . This is 2 to 3 orders of 
magnitudes larger than the experimentally observed N c . 
However, in the preceding analysis we have neglected all 
prefactors, and we emphasize that the result is strongly 
(to the third power) dependent on the experimental pa- 
rameters included in D' ', i.e. humidity, liquid density, 
diffusivity of vapor, and saturated vapor concentration. 



To verify whether the final packing fraction indeed de- 
pends on the spacing between the particles the moment 
the cross-over time is reached, we go back to our ex- 
perimental data. We define the time-dependent pack- 
ing fraction as N(R P /R(t)) 3 . As the droplet evaporates, 
this packing fraction will increase until it reaches its fi- 
nal value $. At the cross-over, defined by A = 1, the 
droplets will have a packing fraction $*. If this <£* is low, 
the amount of crystalline clusters is still very small. On 
the other hand, if $* is high, we expect crystalline clus- 
ters to have formed already. After the cross-over time, 
the interface moves too fast to allow for further ordering, 
and it just presses the ordered particle clusters closer to- 
gether. In figure 5 we show that droplets with a high <£* 
have a high <£: when <£* ^ 0.1 we obtain a final packing 
fraction above the RCP limit. 




will bend radially no more than a particle size. Then, it 
follows by simple trigonometry that the size of a patch 
S will be related to the ball radius Rbaii and the parti- 
cle size R p via: S/Rbaii — arccos(l — R p /Rbaii)- This 
expression predicts a typical patch size of ~ 15/im for a 
ball with Rbaii — 100/ira and R v = 1/im, which is in the 
right order of magnitude as one can observe in figure 2d. 

In conclusion, in this Letter we devise a simple tech- 
nique to create spherical colloidal supraballs relying only 
on droplet evaporation over a robust superhydrophobic 
surface. The supraballs show a highly ordered struc- 
ture if the number of particles inside the drop is large 
enough to trigger early particle clustering. The criti- 
cal number of particles required to obtain an ordered 
particle packing inside the balls depends on the param- 
eters driving the droplet evaporation (through D') and 
the diffusivity of the particles. Hence, by controlling the 
humidity and ambient temperature the supraball pack- 
ing fraction and hence size can be controlled. Massive 
fabrication of micro-compact-supraballs could easily be 
achieved by simply spraying a colloidal solution over the 
micro-structure in a controlled atmosphere. By tuning 
the wetting properties of the particles one could also be 
able to generate the well-known colloidosomes [ ] using 
the same proposed technique. 
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FIG. 5. Final packing fraction <£/ versus the packing fraction 
at the cross-over time <£*. Droplets below a certain <E>* have 
a too low packing fraction at the cross-over time to achieve 
final packing fractions above the RCP limit. The particle 
packing can not only be obtained from the value of the packing 
fraction, but it can also directly be observed from the SEM 
images of the surface of the supraballs. 



We cannot predict the critical <£* theoretically. How- 
ever, we can, retrospectively, use the experimental criti- 
cal <!>* to compute N c . Using that Rbaii/Rp ~ 0.17V 1/3 
at the cross-over, we obtain N c ~ 10 T , which is in the 
same order of magnitude as our experimental results; see 
figure 4b. 

The particle packing in the supraballs can not only be 
assessed by measuring the packing fraction, but it can 
also directly be seen in SEM images from the surface of 
the colloidal supraballs, as shown in figure 5. The size of 
the soccer-ball-like crystalline patches on the surface of 
the ordered supraballs depends on the ball size: bigger 
balls will show larger patches due to the reduced curva- 
ture at their surfaces. To explain the size of the crys- 
talline domains, we hypothesize that a crystalline patch 
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